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Abstract 
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£ — ■ We propose a method to introduce Fock term contributions in relativis- 

tic models of fermions coupled to mesons, including self-interactions for the 
mesonic fields. We show that effects on equilibrium properties and on the dy- 
namical response of the fermionic system can be consistently accounted for. 
| Some implications on equilibrium properties of asymmetric nuclear matter 

are discussed. In particular an indication is emerging for a reduced contribu- 
tion of charged mesons to the symmetry term of the nuclear equation of state 
around normal density. 
P : PACS: 24.10.Cn, 24.10.Jv, 21.30.Fe 
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g I. INTRODUCTION 

Since the pioneering work of Walecka [1[] , the relativistic field model of hadrons, called 
Quantum Hadrodynamics ( QHD ), with its developments and extensions, has been widely 
used in investigations of nuclear systems, see the extensive reviews 0-0. This model can 
be regarded as an effective field theory 0,||,|7j at the nuclear scale, and it is based on a 
relativistic effective Lagrangian which includes strongly interacting nucleon and meson fields. 
The success of QHD has been due to the fact that it provides a realistic description of bulk 
properties both of infinite nuclear matter and of finite nuclei 0-0], even with a minimal 
number of meson fields, the neutral a and u fields and the charged p field. 

The usual procedure to deal with theories which contain baryons interacting with mesons 
consists of eliminating the meson fields in favour of the fermionic degrees of freedom. This 
can be done by exploiting the equations of motion for the meson fields. In this way the 
physical system can be ultimately described in terms of nucleons interacting through the 
exchange of virtual mesons. 

It is important to notice that, once the Fock exchange terms of the nucleon-nucleon 
interaction are taken into account, the Walecka model can reproduce some basic properties 
of nuclear matter by means of the exchange of the scalar a and vector uj fields alone ( see 
Refs. quoted in ). Including the Fock terms also the whole variety of processes in 

nuclear dynamics arising from the fermionic intrinsic degrees of freedom ( spin and isospin ) , 
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can be accounted for |8|- |T0|1 - In principle, if the model is constructed self-consistently, the 
introduction of further meson fields with the corresponding constants for the coupling to 
nucleons appears to be less necessary. We will elaborate in detail on this point. 

However, in its original version, the QHD model suffered from a severe deficiency, namely 
it was not able to give a reasonable value for the nuclear compressibility modulus and the 
effective mass. This unpleasant drawback has been cured by introducing self-interactions of 
the scalar a field up to fourth order |TT| , [T2"| ]. Moreover, within both the modern scheme of 
effective field theories [p!3| , p^| and hadronic chiral models |15| a non-linear self-interaction 
for the scalar meson emerges in a natural way. 

Calculations based on non-linear models of QHD have been so far performed in the 
relativistic mean field approximation ( RMF ), which essentially corresponds to the Hartree 
approximation. With the inclusion of these higher order terms in the Lagrangian, the 
evaluation of the exchange terms of the nucleon-nucleon interaction represents a very difficult 
task. This is due to the non-linearity of the resulting equation of motion for the meson 
field. Although the RMF approximation to QHD with non-linear terms has achieved a large 
success in reproducing quantitatively many observables, in particular for finite nuclei |5|], the 
inclusion of the Fock terms would be highly desirable. In fact, these terms originate from 
the correlations due to the Fermi-Dirac statistics, therefore they are related to a genuine 
quantum effect, which in general cannot be neglected in studying a many-body system. 

In this paper, we present a procedure which allows to evaluate the contribution of the 
Fock terms in a truncation scheme. The result is an expansion of the Fock terms in a 1/N 
series, where N is the multiplicity of the fermionic intrinsic degrees of freedom ( N = 4 for 
nucleons ). In this sense the approach suggested here is not of perturbative nature in the 
coupling constants, which are always accounted for at all orders. 

We illustrate our method in the case of infinite nuclear matter. We show that both 
equilibrium properties and dynamics can be self-consistently treated at the same degree of 
approximation. 

We would like to remark that our treatment of the Fock terms is not restricted to QHD 
models. It can be applied in general, most likely with non trivial consequences, to models 
involving fermions coupled to bosons. Among these, we mention the Quark Meson Coupling 
model where nuclear matter is described as a collection of composite nucleons with 



meson fields directly coupled to quarks |fPlj-|T!|. Outside a strictly nuclear context, we 



mention also the Linear Quark Meson model, which represents an effective description of 
QCD below the mesonic compositeness scale ( ~ 600MeV ) ( see the reviews []2U| and 
This model provides a hybrid description of QCD in terms of quarks and mesons. 



II. TREATMENT OF FOCK TERMS IN A KINETIC APPROACH 

The aim of the present paper is to illustrate a procedure which allows to take into account 
the self-interactions of meson fields within a Hartree-Fock scheme. For simplicity, and also 
to focus on the new procedure, we limit ourselves to consider only the minimal number of 
mesons which is sufficient to reproduce some basic properties of nuclear matter. Therefore 
our starting point is the QHD-I a — uj model D, including self-interaction terms for the 



scalar a field 11,12|. The Lagrangian density for this model is given by: 
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C = $bv{id» - g v Vn — (M — g s m + 5 ( d ^<t> ~ m s 2 <t> 2 ) 

-f^ ~ h ~/ - \W, V W^ + l -m v ^ , (1) 

where W JJ/ (x) = d^V v {x) — d v V^{x) . Here ip(x) denotes the 8(spin and isospin)-component 
nucleon field, the scalar (<fi(x)) and the vector (V u (x)) boson fields are associated with the 
scalar (a) and vector ( u ) mesons, respectively (h — c — 1). 

Fock contributions naturally appear in the source terms of the equation of motion of 
the fields. Here we will study the related effects on the corresponding kinetic equations, in 
particular for the fermion field. This will allows us to derive the contributions of exchange 
terms on static and dynamical properties of the interacting nuclear system in a self-consistent 
way, at the same level of approximation. 



A. The Wigner transform formalism 

In agreement with the Relativistic Mean Field (RMF) picture of the QHD model we 
focus our analysis on a description of the many-body nuclear system in terms of one-body 
dynamics. Therefore, it is convenient to use the Wigner transform of the one-body density 
matrix of the fermionic field. The one-particle Wigner function is defined as: 



where a and (3 are double indices for spin and isospin. The brackets denote statistical 
averaging and the colons denote normal ordering. According to the Clifford algebra, the 
Wigner function F(x,p) can be decomposed as: 

F(x,p) = F(x,p) + j^(x,p) + ^a^F^(x,p) 

+ j 5 A(x,p) -7 5 7^ M (x,p) . (2) 
The integral over p of the Wigner function 

[F{x)\ a i3 = / d A p [F(x,p)] a p = <: $p{x)$ a {x) :> 



is related to the various densities with their specific transformation properties both in or- 
dinary space and in isospin space ( isoscalar and isovector ). For instance the scalar and 
current isoscalar densities are given by: 

Ps(x) = <: ip(x)ip(x) :>= TrF(x) = 8 F(x) , 

f(x) = <: tp(x)^ij(x) :>= TrYF(x) = 8 F fl (x) . 

The equation of motion for the Wigner function can be derived from the Dirac equation 
by using standard procedures (see e.g. Refs. |2"2|]2"3]]): 
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9V77^TI f d 4 Re'^ R <: Mx + )^ a5 M^-W,(x-) :> 



+9st^ti I d 4 Re~ ip - R <: ^(ar + )^«(a:_)0(x_) :>= , (3) 



(2vr) 4 

1 

with x + = x + § and x_ = x — 

B. The truncation scheme for exchange contributions 

The statistical averages of the fermionic and mesonic operators in the above equation 
will now be treated within a Hartree-Fock ( HF ) approximation. In order to take into 
account the contribution of exchange terms in a manageable way we assume, as a basic 
approximation, that in the equations of motion for the meson fields the terms containing 
derivatives can be neglected with respect to the mass terms. Then, the meson field operators 
are connected to the operators of the nucleon scalar and current densities by the following 
equations: 

$/f s + + = ^(x)^(x) = ps(x) , 

V»(x)=f v 4,(x)rHx)=Ux), (4) 

where $ = g s <f), V" = g v V» and f s = (g s /m s ) 2 , f v = (g v /m v ) 2 , A = a/g% B = bjg%. 

This approximation is not justified for light mesons, like pions. The inclusion of self- 
interaction terms of the pionic field is a very difficult task and a different approximation 
scheme is needed. However in this case a perturbative expansion in the pion-nucleon coupling 
constant seems to be reasonable [^|. Moreover it has been shown that the inclusion of pions 



does not change qualitatively the description of nuclear matter around normal conditions 



24 



The next step will be the elimination of the meson fields in Eq.(3) by means of Eq.(4). 
As far as the term containing the vector field is concerned, by a simple substitution one 
can obtain a statistical average of only fermionic field operators ||. Instead, because of 
the non-linear coupling of the a field to the fermionic scalar density, the elimination of the 
scalar field in the last term of Eq.(3) requires a particular procedure. 

We consider the expansion of the field operator $ in terms of the scalar density operator 
Ps'- & — J2nXnPs- Then, in the HF approximation , by means of the Wick's theorem the 
average in the last term of Eq. (3) can be expressed as: 

<: if}{3(x + )il} a (x-)§(x-) :>=<: ^p(x+)if} a (x-) :><: :> + <: ^ /3 (x + )^ 7 (a;_) :> 

The statistical averages of and its derivatives are given by series of <: Ps(x) :>, which 
in turn can be expanded as: 
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<: :>= pS(x) - £ ^"^^T^ " (6) 

As a consequence, we obtain for <: :> the expansion: 

<=*(,) :>=«(,)-gI(-l)'7V^(,)gg, (7, 

where = Y,nXnPs( x )i an d similar expressions for <: d k §(x)/d k ps(x) :>. The expan- 

sions of Eqs.(5),(6) and (7) allow to settle a truncation scheme. The parameter to fix the 
approximation order is given by the number of factors ps{x) = ip(x)i(j(x) which are broken 
in writing Eq.(5) with the expansion (7). For each decoupling there is a trace operation 
less, leading each time to a quenching factor of 1/8. Moreover in our case the coefficients A 
and B in Eq.(4), to be determined by fitting the bulk properties of nuclear matter, should 
be small |0| . Therefore, the derivatives d k Q(x) / dp k s (x) become smaller and smaller as k 
increases. 

The expansion Eq.(5) actually implies that contributions of nucleons of the Dirac sea 
are neglected. This approximation will allow us to derive a closed equation for the Wigner 
function containing only finite quantities. We discuss this point at the end of the section. 

We consider exchange corrections only up to the next-to-leading term, with respect to 
the Hartree contribution. In this approximation Eq.(5) assumes the form: 

- 1 d 2 <fr{x) 

<: i)f3(x + )ip a (x^)^(x-) :>=<: ij)p{x + )il) a (x-) :> - -TrF 2 (x. N 



dPs( x -) 

- <: ^(x + )^ 7 (z_) :> [F(x-)] *®y~\ , (8) 

dp s {x-) 

where the function obeys the equation: <&(x)/fs + A& 2 (x) + 5$ 3 (x) = ps(x). 

According to the procedures usually followed in the derivation of kinetic equations (Refs. 
BHH) (semi-classical approximation), we assume that the densities <: ip^(x^)ip a (x^) :> and 
the field $(a;_) are slowly varying functions, and we retain only the first order term in their 
Taylor expansion at the point x. This is essentially the same kind of assumption previously 
done for the meson fields and expressed by Eq.(4). In fact because of the small Compton 
wave-lengths of the heavy mesons a and lu, retardation and finite range effects can be 
neglected. 

Within the approximations outlined above Eq.(3) becomes 

Z -d^F(x,p) + Y(P» ~ fvUx))F(x,p) -(M- $(x) + ^TrF 2 (x)^^jF(x,p) 



■5 A ^)r-*(*) + i^)gg 



M^A - l)YF(x) lfl F(x,p)] + [(~A - = > ( 9 ) 



where A = d x ■ d p , with d x acting only on the first term of the products. 
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This is the kinetic equation for the matrix F(x,p) in mean-field dynamics. In addition 
to the scalar and vector isoscalar fields, the exchange contributions, which are the terms 
containing the matrix F(x) in Eq.(9), may give rise to scalar, vector, tensor, pseudoscalar 
and pseudovector fields ( both isoscalar and isovector ), leading to completely new effects. 

In Eq.(9) we have a general expression for the effective mass, including an isospin contri- 
bution. Since the derivatives of mesonic fields have been neglected, the effective mass does 
not present a momentum dependence. 

A quantity of interest in the study of nuclear dynamics is the statistical average of the 
canonical energy-momentum density tensor. Since terms containing the derivatives of the 
meson fields are neglected in our approximation, the tensor takes the form: 

= l ^{x) llx d^{x) + [U@) - ^-V x (x)V\x)]g^ (10) 

where g^ v is the diagonal metric tensor and U(<§>) = / fs + ^4/3 $ 3 + B/A $ 4 . According 
to the approximation introduced above, the energy-momentum density tensor, in the HF 
approximation, reads: 

<: T^(x) :>= 8 J d 4 p p v F^x,p) + {[/($) - ^j x ( x )j x (x)}g^ 

~l [Trp2{x) l^ " /vTr( 7A F(x) 7 A F(x))]^ (11) 

The quantities in square brackets are the exchange contributions. It is essential to note that 
Fock terms contain traces of powers of F(x) that naturally may bring contributions from all 
the fields having tensorial properties consistent with the symmetry of the system. 

/,From Eq.(ll) we obtain the energy density for symmetric nuclear matter that in analogy 
to the Hartree case can be rewritten in the following form: 

e = <: Too :> = J dpE;9(p F - |p|) + U($) + ^f sPs 2 + \f V p B \ (12) 



where E* = y |p| 2 + M* 2 , ps = jo is the baryon density and p F is the Fermi momentum. 
M* is the nucleon effective mass: 



1 . 1 d$ , 1, 2 2 * d 2 $ 



M' = M-*-(- !v ---) Ps + T y s+A)d 



The scalar density ps is given by 



The density dependent coupling constants fs and f v are defined as it follows: 

fs -2 fv -% { dp- s +Ps Jpl )] 
7 5 , l.d$ d 2 <I> 

f v = 4 fv + -^w s - ps M ] - (13) 
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To give some indications about the convergency of the truncation in Eq.(8), we plot in 
Fig.l the density behaviour of the coupling functions fs and fv, obtained using for f s , f v , 
A and B the parameterization given in chapter III. The dashed lines represent the results 
obtained neglecting the terms containing the second derivative of $ in Eq.(13), while the 
full lines correspond to the complete calculations (Eq.(13)). Up to densities of the order 
of 4 po, it is possible to observe that, for the parameterization considered, the inclusion of 
the second derivatives changes slightly the results. Therefore in such a density domain one 
expects a even smaller contribution from the next terms of the expansion of Eq.(6), since 
they contain higher order derivatives. 

In order to establish a link with previous works on the introduction of Fock terms in 
QHD theories, we will consider only the linear term in the equation for the scalar field, 
Eq.(4), where this kind of calculations have been performed so far [10,24-27], [f28|j . 

Our results are clearly reducing to the ones reported in Ref.[8]. The single-particle 
energies and the nucleon effective mass are given by: 

e p = e; + f vPB + [(i/ v + lf 8 ) PB ] = e; + f vPB (w) 



and 



M* = M- f s ps - i(\fv - \fs)ps] =M-(f s + f s )p s , (15) 

respectively, where now fs and fv of Eq.(13) are simply given by: fs = —gfs + \ fv , 
fv = hfs + ifv ■ The energy density, e =<: T 00 :>, is expressed by: 



(27T) 



J d P E;e(p F - \ P \) + ±fsp s + \fy P \ 



+ \[{\fv-\fs)p 2 s + {\fv + \fs)p 2 B } 



(27T) 



J dpE* p 9(p F - |p|) + \{fs + ~fs)Ps + \fvpl, (16) 



with the exchange contributions in Eqs.(14),(15) and (16) given by the terms in square 
brackets. 

We remark that the same expressions for E p , M* and e can be obtained from the results 
of the HF calculations of Ref. P^J neglecting the square of the four-momentum transfer 



between two interacting nucleons with respect to the squared meson masses. Hence for 
linear models of QHD our approach is equivalent to the usual HF approximation, once 
retardation and finite-range effects are neglected. 

The two-loop approximation for nuclear matter was introduced to study the loop 
expansion scheme for QHD. According to this approximation the contribution of exchange 
diagrams to the energy density is evaluated using nucleon propagators in the Hartree approx- 
imation (i.e. at one loop level). Neglecting vacuum, retardation and finite range effects, the 
results of Ref. [25] can be reproduced by using a similar approximation in our approach. This 
amounts to neglect the exchange terms in the equation for the equilibrium Wigner function, 
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and to retain the exchange contributions in the equation for the energy density, Eq.(16). 
Therefore, within this approximation, the exchange terms in the nucleon self-energy are 
disregarded. However it should be noticed that their contribution to the energy density, i.e. 
the second term in the square bracket of Eq. (16), Ae = (\fv + \f s )p 2 Bi is comparable to 
the exchange contributions that are included in |25] and hence it should not be neglected . 
A similar situation occurs in the self-consistency equation for the nucleon effective mass. 

Within the formalism used in the present paper, vacuum effects could be taken into 
account by omitting the normal ordering and subtracting a vacuum term in the definition 
of the Wigner function |23|| . Such a subtraction does not eliminate all the divergent vacuum 
terms in the equations for the Wigner function and the energy-momentum density tensor, 
thus a renormalization procedure in necessary [23||. Since, as stressed before, our approach 
is equivalent to the HF approximation, we should expect to get results analogous to those 
obtained in |25|j27|1 . A fully self-consistent HF description of QHD including vacuum effects 
is very complicated |27j and has not been achieved yet. On the other hand, two-loop 
corrections to the Hartree approximation give unnaturally large vacuum contributions to 
the energy of nuclear matter |25| . This would suggest to follow other paths for representing 
the vacuum dynamics. For a review on these points see also Ref. [Q. 

However, it is generally recognized that the RMF approximation, where negative energy 
states of nucleons are not considered, is very successful in describing properties both of 
nuclear matter and of finite nuclei, once self-interaction terms of the a-field are included. 
The approach presented here can be considered as an extension of the RMF approximation 
to include exchange terms on the same basis. 



III. INFLUENCE OF FOCK TERMS ON THE NUCLEAR EQUATION OF STATE 

The study of the nuclear Equation of State (EOS) can be performed by evaluating the 
statistical average of the energy-momentum tensor, Eq. (11). 

Specializing to the case of symmetric nuclear matter, we have tested that our procedure 
to introduce the Fock terms leads to a thermodynamically consistent theory. The equilibrium 
Wigner function together with the single-particle energy spectrum can be obtained following 
the scheme of Ref. ||. We have found that the relation between pressure and energy density: 

p ="^- t = H- :> (17) 



and the Hugenholtz-Van Hove (HV) theorem [29 



de 

H= — = e F (18) 



dp 



B 



are satisfied. Here e =<: T 00 :>, pb — j° is the baryon density and ep is the Fermi energy. 
This actually represents a good check of the consistency in the truncation procedure used. 
We remind that the explicit inclusion of meson field derivatives will make more delicate the 
proof of the HV theorem |30] . 

In order to illustrate how the properties arising from the nucleonic intrinsic degrees of 
freedom can be described with our approach, we now present some results about the EOS 
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of asymmetric nuclear matter at equilibrium. The matrix F a p{x,p) can be decomposed in 
the isospin space as: F a/3 (x,p) = {F s ) a ipi{x,p) + T3(F 3 ) ,i^i(x,p), where r 3 is the third Pauli 
matrix, (a/3) are spin-isospin indices, while {u'f3') denote spin indices. In symmetric nuclear 
matter the matrix F 3 (x,p) vanishes, whereas in the asymmetric case its components give the 
various isovector densities. Due to the asymmetry of the system, the kinetic energy density 
increases because the Fermi momenta of protons and neutrons are different. The presence 
of the exchange terms (terms in square bracket in Eq.(ll)) gives rise also to a contribution 
from the potential part of the energy density e. 

1 <9 2 e, 

The symmetry energy can be defined as E svm = — — — \ I=0 , where / is the asymmetry 

p B oP 

parameter, defined as the ratio {{pb)u— (pb) p )/pb (n stands for neutrons and p for protons). 
In Fig.l we have reported the calculated symmetry energy as a function of the baryon density 
Pb- In the calculations the four parameters of the model (fs, fv, A and B) have been fitted 
in order to reproduce equilibrium properties of symmetric nuclear matter: saturation density 
Po = 0.16/m -3 , binding energy E/A=-15.75 MeV, nucleon effective (or Dirac) mass at p , 
Mq* = 0.7M, incompressibility modulus K = 250MeV. In this way we obtain the following 
value for the parameters: fs = 10.1, fy = 4.1, A = 0.06 and B = 0.01. In the scalar field 
contribution we clearly see a steady decrease of the weigth of higher order terms. This 
seems to be a qualitative signal of a naturalness in the expansion. However, a more accurate 
discussion should be done only in a model that includes non-linearly both scalar and vector 
fields p|,|31 



A quite interesting feature of our approach is that now the asymmetric nuclear matter 
is described relatively well with the same parameters fixed by the fit. Figure 2 shows that 
our non-linear HF (NLHF) evaluation (solid line) gives a symmetry energy value at normal 
density, i.e. the coefficient of the Weiszaecker mass formula of = 24 MeV. This result 
appears too low with respect to the accepted value (around 28 — 30MeV), but we have to 
consider that only the contribution of isoscalar mesons are included. From this point of 
view our calculation shows the importance of exchange terms on the symmetry energy, even 
though the inclusion of the p meson is necessary to achieve a good value of 04. 

To be more quantitative, considering that the kinetic contribution to the (Z4 parameter 
is of the order of YIMeV (~ c,f/3), of the remaining part of about VoMeV almost two 
thirds seem to come from exchange terms. This appears to be an interesting indication of a 
reduced contribution of the p meson strength in nuclear matter around normal density. 

A rather repulsive density dependence is also deduced (Fig.2; full line). Non-linear 
Hartree (NLH) calculations including the p meson (dashed line) can reproduce the same 04 
value once the coupling constant f p is suitably fitted. However in the NLHF calculations 
we also obtain a non-linear density dependence of the symmetry energy that arises from 
the non-linear scalar self-interactions through exchange terms. This effect is not present in 
NLH calculations. 



IV. CONCLUSIONS 

We have introduced a procedure to investigate the role of the Fock exchange terms in a 
relativistic model of fermions coupled to mesons, including self-interactions for the mesonic 
fields. We have shown that the evaluation of the Fock terms, within a suitable truncation 
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scheme, leads to a consistent description of the equilibrium and dynamical properties related 
to the fermionic intrinsic degrees of freedom. An application to the case of asymmetric 
nuclear matter at equilibrium gives satisfactory results concerning the density dependence 
of the symmetry energy. It should be noticed that also other properties, such as the isospin 
and density dependence of effective masses, that can be deduced from the stationary solution 
of the kinetic equation, are reasonably reproduced within our procedure, in comparison to 
QHD models including also scalar isovector mesons P^fl . A quite interesting general feature 
is appearing: a consistent inclusion of exchange contributions seems to clearly reduce the 
weight of charged mesons in describing isospin properties of nuclear matter around and 
above normal barion density. 

Finally we remark that the aim of this work is mainly to present a procedure to evaluate 
Fock terms of a non-linear field coupling scheme. For this reason we have omitted a full 
analysis of all possible features of a EOS calculation; in particular our model is subject 
to the possible criticism of all the approaches with non-linearities only in the scalar field. 
However it is straightforward to extend the same procedure to non-linear vector fields. Thus 
we stress that the method described here is very general and can be applied to effective field 
theories not only for nuclear matter but in general to systems with fermion-boson dynamical 
couplings. Of course, the analysis can be extended also to the investigation of spin effects. 
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Figure Caption 



Fig.l: Density behaviour of the coupling functions fs (bottom panel) and /y (top panel). 
See the text for details. 

Fig. 2: Symmetry energy per nucleon vs. baryon density. Full line: present NLHF 
calculation. Dashed line: NLH calculation including the p meson. 
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